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Abstract. High-accuracy astrometry permits the determina- 
tion of not only stellar tangential motion, but also the compo- 
nent along the line-of-sight. Such non-spectroscopic (i.e. astro- 
metric) radial velocities are independent of stellar atmospheric 
dynamics, spectral complexity and variability, as well as of 
gravitational redshift. Three methods are analysed: (1) chang- 
ing annual parallax, (2) changing proper motion and (3) chang- 
ing angular extent of a moving group of stars. All three have 
significant potential in planned astrometric projects. Current 
accuracies are still inadequate for the first method, while the 
second is marginally feasible and is here applied to 16 stars. 
The third method reaches high accuracy (< 1 km s _1 ) already 
with present data, although for some clusters an accuracy limit 
is set by uncertainties in the cluster expansion rate. 

Key words: Methods: data analysis — Techniques: radial ve- 
locities — Astrometry — Stars: distances — Stars: kinematics 
— open clusters and associations: general 



1. Introduction 

For well over a century, radial velocities for objects outside 
the solar system have been determined through spectroscopy, 
using the (Doppler) shifts of stellar spectral lines. The advent 
of high-accuracy (sub-milliarcsec) astrometric measurements, 
both on ground and in space, now permits radial velocities to 
be obtained by alternative methods, based on geometric prin- 
ciples and therefore independent of spectroscopy. The impor- 
tance of such astrometric radial velocities stems from the fact 
that they are independent of phenomena which affect the spec- 
troscopic method, such as line asymmetries and shifts caused 
by atmospheric pulsation, surface convection, stellar rotation, 
stellar winds, isotopic composition, pressure, and gravitational 
potential. Conversely, the differences between spectroscopic 
and astrometric radial velocities may provide information on 
these phenomena that cannot be obtained by other methods. 
Although the theoretical possibility of deducing astrometric ra- 
dial velocities from geometric projection effects was noted al- 
ready at the beginning of the 20th century (if not earlier), it is 



only recently that such methods have reached an accuracy level 
permitting non-trivial comparison with spectroscopic measure- 
ments. 

We have analysed three methods by which astrometric ra- 
dial velocities can be determined (Fig. [[]). Two of them are ap- 
plicable to individual, nearby stars and are based on the well 
understood secular changes in the stellar trigonometric parallax 
and proper motion. The third method uses the apparent changes 
in the geometry of a star cluster or association to derive its kine- 
matic parameters, assuming that the member stars share, in the 
mean, a common space velocity. In Sects. ^ to ^ we describe 
the principle and underlying assumptions of each of the three 
methods and derive approximate formulae for the expected ac- 
curacy of resulting astrometric radial velocities. For the first 
and second methods, an inventory of nearby potential target 
stars is made, and the second method is applied to several of 
these. 

However, given currently available astrometric data, only 
the third (moving-cluster) method is capable of yielding as- 
trophysically interesting, sub-km s _1 accuracy. In subsequent 
papers we develop in detail the theory of this method, based 
on the maximum-likelihood principle, as well as its practical 
implementation, and apply it to a number of nearby open clus- 
ters and associations, using data from the Hipparcos astrometry 
satellite. 



2. Notations 

In the following sections, n, /.i and v r denote the trigonomet- 
ric parallax of a star, its (total) proper motion, and its radial 
velocity. The components of /i in right ascension and declina- 
tion are denoted and fig, with /i = (//^ + Hg) 1 ^ 2 - The 
dot signifies a time derivative, as in tt = dir/dt. The statisti- 
cal uncertainty (standard error) of a quantity x is denoted e(x). 
(We prefer this non-standard notation to e x , since x is itself 
often a subscripted variable.) a v is used for the physical ve- 
locity dispersion in a cluster. A = 1.49598 x 10 8 km is the 
astronomical unit; the equivalent values 4.74047 km yr s _1 
and 9.77792 x 10 s mas km yr s _1 are conveni ently used in 
equations below (cf. Table 1.2.2 in Vol. 1 of ESA |l997| ). Other 
notations are explained as they are introduced. 



* Based (in part) on observations by the ESA Hipparcos satellite 
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3. Astrometric accuracies 

In estimating the potential accuracy of the different methods, 
we consider three hypothetical situations: 

- Case A: a quasi-continuous series of observations over a 
few years, resulting in an accuracy of e(ir) = 1 mas 
(milliarcsec) for the trigonometric parallaxes and e(/i) = 
1 mas yr -1 for the proper motions. 

- Case B: similar to Case A, only a thousand times better, i.e. 
e(7r) = 1 /ias (microarcsec) and e(/i) = 1 fias yr -1 . 

- Case C: two sets of measurements, separated by an inter- 
val of 50 yr, where each set has the same accuracy as in 
Case B. The much longer-time baseline obviously allows a 
much improved determination of the accumulated changes 
in parallax and proper motion. 

The accuracies assumed in Case A are close to what the 



Hipparcos space astrometry mission (ESA 1997) achieved for 
its main observation programme of more than 100000 stars. 
Current ground-based proper motions may be slightly better 
than this, but not by a large factor. This case therefore repre- 
sents, more or less, the state-of-the-art accuracy in optical as- 
trometry. Accuracies in the 1 to 10 fias range are envisaged 
for some planned or projected space astrometry missions, such 



as GAIA (Lindegren & Perryman 1996) and SIM (Unwin et 
al. 1998| ). The duration of such a mission is here assumed to 
be about 5 years. Using the longer-time baselines available 
with ground-based techniques, similar performance may in the 
future be reached with the most accurate ground-based tech- 
niques (Pravdo & Shaklan |1996| ; Shao |1996j ). Case B there- 
fore corresponds to what we could realistically hope for within 
one or two decades. Case C, finally, probably represents an 
upper limit to what is practically feasible in terms of long- 
term proper-motion accuracy, not to mention the patience of 
astronomers. 



4. Radial velocity from changing annual parallax 

The most direct and model-independent way to determine ra- 
dial velocity by astrometry is to measure the secular change 
in the trigonometric parallax (Fig. [j]a). The distance b (from 
the solar system barycentre) is related to parallax tt through 
b = A/ simr ~ A/ir. Since v r — b, the radial velocity is 



-A- 



(1) 



where A is the astronomical unit (Sec t. pf). The equivalent of 
Eq. (Tij) was derived by Schlesinger (1917), who concluded 



that the parallax change is very small for every known star. 
However, although extremely accurate parallax measurements 
are obviously required, the method is not as unrealistic as it 
may seem at first. To take a specific, if extreme, example: for 
Barnard's star (Gl 699 = HIP 87937), with tt = 549 mas 
and v r — —110 km s -1 , the expected parallax rate is tt = 
+34 /ias yr -1 . According to our discussion in Sect. |^ this will 
almost certainly be measurable in the near future. 



(a) Changing annual parallax 
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Fig. 1. Three methods to astrometrically determine stellar ra- 
dial motion: a Changing trigonometric parallax tt; b Perspec- 
tive change in the proper motion /i; c Changing angular separa- 
tion p of stars sharing the same space velocity, e.g. in a moving 
cluster. 



It can be noted that the changing-parallax method, in con- 
trast to the methods described in Sect. |]and |J does not depend 
on the object having a large and uniform space motion, and 
would therefore be applicable to all stars within a few parsecs 
of the Sun. 

4.1. Achievable accuracy 

The accuracy in v r is readily estimated from Eq. ([j]) for a given 
accuracy in tt, since the contribution of the parallax uncertainty 
to the factor A/tt 2 is negligible. The achievable accuracy in tt 
depends both on the individual astrometric measurements and 
on their number and distribution in time. Concerning the tem- 
poral distribution of the measurements we consider two limit- 
ing situations: 
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Table 1. Number of target stars, and achievable accuracy in 
the radial velocity, as obtained from the changing annual par- 
allax, n = number of stars with parallax greater than it in 



the Third Catalogue of Nearby Stars (Gliese & JahreiB |1991 ); 
e(v r ) = predicted accuracy according to Eqs. (fy and (|3). 
Case B: space astrometry mission lasting L = 5 yr and yielding 
parallaxes with standard error e(-7r) — 1 /ias and consequently 
e(7r) = 0.69 /ias yr -1 ; Case C: combination of two such mea- 
surements, e(7Ti) = e(7r 2 ) = 1 /ias, with a large epoch differ- 
ence, T = 50 yr. 
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Quasi-continuous observation. The measurements are 
more or less uniformly spread out over a time period of length 
L centred on the epoch to. This is a good approximation to 
the way a single space mission would typically be operated; 
for example, Hipparcos had L ~ 3 yr and to ~ 1991.25. In 
such a case there exist simple (mean) relations between how 
accurately the different astrometric parameters of the same 
star can be derived, depending on L. For instance, e(fig) — 
(Vl2/ L)e(8), if e{ns) is the accuracy of the proper motion in 
declination and e(<5) that of the declination at to- This approxi- 
mation is applicable to Case A and B as defined in Sect. 

Two-epoch observation. Two isolated parallax or proper- 
motion measurements are taken, separated by a long time in- 
terval (say, T years) during which no observation takes place. 
Each measurement must actually be the result of a series cov- 
ering at least a year or so, but the duration of each such series 
is assumed to be negligible compared with T. This could be 
two similar space missions separated by several decades and is 
applicable to Case C in Sect. [|. 

For quasi-continuous observation we may assume that the 
parallax variation is linear over the observation period L. Thus, 
ir(t) = ttq + (t — to)7r, where tto and 7r are two parameters to be 
determined from the observations. There exists an approximate 
relation between the accuracies of these two parameters that is 
similar to that between the proper motion and the position at 
the mean epoch, viz. e(n) — (Vl2/ L)e(no). Moreover, the es- 
timates of the two parameters are uncorrelated, so c(tto) equals 
the accuracy e(7r) of a parallax determination in the absence of 
the parallax-change term; thus 



e(v r ) 



12A 



Ltt 2 



(2) 



In the case of a two-epoch observation, let us assume that inde- 
pendent parallax measurements tti and 7r 2 are made at epochs 
ti and < 2 = *i + T. The estimated rate of change is tt = 
(7r 2 — 7r 1 )/T. With e(7r 1 ) and e(7r 2 ) denoting the accuracies of 



the two measurements, we have t{-k) = [e(7r 1 ) 2 +e(7r 2 ) 2 ] 1 / 2 /T 
and consequently 



e(v r ) 



A 
TV 2 " 



[e(7n) 2 +e(7r 2 ) 



21 1/2 



(3) 



For given observational errors we find, from both Eq. (||) and 
(|^), that the radial-velocity error is simply a function of dis- 
tance. The number of potential target stars for a certain maxi- 
mum radial-velocity uncertainty is therefore given by the total 
number of stars within the corresponding maximum distance. 
Table [l] gives the actual numbers of such stars, and the obser- 
vational accuracies that may be reached. 

5. Radial velocity from changing proper motion 
(perspective acceleration) 

To a good approximation, single stars move with uniform lin- 
ear velocity through space. For a given linear tangential veloc- 
ity, the angular velocity (or proper motion /i), as seen from the 
Sun, varies inversely with the distance to the object. However, 
the tangential velocity changes due to the varying angle be- 
tween the line of sight and the space- velocity vector (Fig. |l]b). 
As is well known (e.g. van de Kamp 1967, Murray 1983| ) the 
two effects combine to produce an apparent (perspective) ac- 
celeration of the motion on the sky, or a rate of change in proper 
motion amounting to /i = — 2/iu,./fr. With b = A/ir we find 



-A- 



27T/1 



(4) 



Schlesinger ( |1917[ ) derived the equivalent of this equation, 
calculated the perspective acceleration for Kapteyn's and 
Barnard's stars (cf. Table ||) and noted that, if accurate posi- 
tions are acquired over long periods of time, "we shall be in 
position to determine the radial velocities of these stars inde- 
pendently of the spectroscope and with an excellent degree 
of precision". The equation fo r the perspective acceleration 
was ea rlier derived by Seeliger ( 190l| )F| and used by Ristenpart 
( 1902 ) in an (unsuccessful) attempt to determine fi observa- 
tionally for Groombridge 1830. A major consideration for Ris- 
tenpart seems to have been the possibility to derive the paral- 
lax from the apparent acceleration in combination with a spec- 
troscopic radial velocity. Such a determination of 'acceleration 
parallaxes' was also considered by Eichhorn (1981). 

Subsequent attempts to determine the perspective accelera- 
tion of Barnard's star by Lundmark & Luyten (1922), Alden 
(1924) and van de Kamp (1935b) yielded results that were 
only barely significant or (in retrospect) spurious. Meanwhile, 
Russell & Atkinson (1931 ) suggested that the white dwarf van 



Some remarks in the literature, e.g. by Ristenpart (1902) and 



Lundmark & Luyten (1922), seem to s ugges t that the perspective ac- 
celeration was discovered by Bessel (1844). However, as far as we 



can determine, Bessel only discussed proper-motion changes caused 
by gravitational perturbations, explicitly neglecting terms depending 
on the radial motion. 
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Table 2. Nearby high-proper-motion stars suitable for the determination of radial velocity from changing proper motion (per- 
spective acceleration). All known objects with a parallax-proper-motion product greater than 0.5 arcsec 2 yr -1 are included. Data 



are from the Hipparcos Catalogue (ESA 1997) where available, otherwise from the preliminary version of the Third Catalogue of 
Nearby Stars, CNS3 (Gliese & JahreiB |1991| ). Columns: CNS3 = identifier in CNS3; HD = identifier in the HD/HDE catalogue; 
HIP = identifier in the Hipparcos Catalogue; Sp = spectral classification in CNS3; V = visual magnitude; b = distance from the 
Sun; 7r/i = product of parallax and proper motion; e(v r ) = predicted accuracy of the astrometric radial velocity (Case B: space 
astrometry mission lasting L = 5 yr and yielding an accuracy of e(/i) = 1 /ias yr -1 for the mean proper motion and consequently 
e(/i) = 1.55 //as yr -2 for the acceleration; Case C: the combination of two such measurements with e(//i) = e(^2) = 1/xasyr -1 
having an epoch difference T = 50 yr). In the last column, P is the orbital period of a binary. Possibly, the star CF UMa does 
not exist, being an erroneous catalogue entry. 
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Barnard's star 
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M5Ve 
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0.25 


0.01 


a Cen C (Proxima) 
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a Cen A (AB: P = 80 yr) 
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Kapteyn's star 
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61 Cyg B (AB: P = 700 yr) 
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UV Cet (AB: P = 27 yr) 
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Luyten's star 
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P = 2.2 yr 
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40 Eri B (BC: P = 250 yr) 


Gl 166C 






dM4.5e 


9.5 


4.8 


0.84 


0.90 


0.03 


DYEri 


G1299 






dM5 


12.8 


6.8 


0.77 


0.98 


0.04 


Ross 619 


G1451A 


103095 


57939 


G8VI 


6.4 


9.2 


0.77 


0.98 


0.04 


Groombridge 1830 


G1451B 








12 


9.2 


0.82 


0.92 


0.03 


CF UMa (non-existent star?) 
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AB: P = 20 yr 



Maanen 2 might exhibit a gravitational redshift of several hun- tion. The astrophysical relevance of astrometric radial-velocity 
dred km s _1 and that this could be distinguished from a real ra- determinations was thus already established (Oort 1932). 
dial velocity through measurement of the perspective accelera- In relatively recent times, the perspective acceleration was 

successfully determined for Barnard's star by van de Kamp 
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(1962, 1963, 1967, 1970, 1981); for van Maanen 2 by van 



deKamp (1971), Gatewood & Russell (1974) and Hershey 



(1978); and for Groombridge 1830 by Beardsley et al. (1974) 



Among these determinations the highest precisions, in terms 
of the astrometric radial velocity, were obtained for Barnard's 
star (corresponding to ±4 km s ~ 1 ; van de Kamp 1981) and van 



Maanen 2 (±15 km s _1 ; Gatewood & Russell 11974 



Our application of the method, combining Hipparcos mea- 
surements with data in the Astrographic Catalogue, yielded ra- 
dial velocities for 16 objects, as listed in Table 

5.7. Achievable accuracy 

The accuracy of the radial velocity calculated from Eq. (Q) can 
be estimated as in Sect. 4.1. It depends on the parallax-proper- 
motion product TT/.1. The most promising targets for this method 
are listed in Table ^, which contains the known nearby stars 
ranked after decreasing tt/i. 

For quasi-continuous observation during a period of length 
L we may use a quadratic model for the angular position <f> of 
the star along the great-circle arc: 4>(t) = 4>(to) + (t — fo)Mo + 
§(t — to) 2 A- Here no is the proper motion at the central epoch 
to - The estimates of no and (i are found to be uncorrected and 
their errors related by e(/i) = (\/60/L)e(/xo). Consequently, 



e(v r ) 



15 A 



LlTfJ, 



(5) 



where e(/i) is the accuracy of proper-motion measurements in 
the absence of temporal changes. We neglect the (small) con- 
tribution to e(v r ) from the uncertainty in the denominator 7r/x. 

For a two-epoch observation, consider proper-motion mea- 
surements n\ an d ^2 made around t\ and t% = ti + T. The 
estimated acceleration is (i — (/i 2 ~ f 1 i)/T- Provided the 
two observation intervals centred on t\ and t% do not overlap, 
the measurements are independent, yielding the standard error 
e(p) = [e{ni) 2 + e(> 2 ) 2 ] 1/2 /T. For the radial velocity this 



gives 



e(v r ) 



A 



\e(ni) 2 + e(/i 2 ) 



21 V2 



(6) 



Based on these formulae, Table || gives the potential radial- 
velocity accuracy for the two cases B and C defined in Sect. 0. 

In a two-epoch observation we normally have, in addition, 
a very good estimate of the mean proper motion between t\ 
and t%, provided the positions <j>\ and </> 2 at these epochs are 
accurately known. In the previous quadratic model we may take 
the reference epoch to be to = (ti + t-x)l2 and find no = {<p2 — 
<j>i)/T with standard error e(po) = [e(0i) 2 + e((f> 2 ) 2 } 1/2 /T. 
The three proper-motion estimates no* Ml an d M2 (referred to 
to, ti and t 2 ) are mutually independent and may be combined 
in a least-squares estimate of /i. If e(ni) = e(/i2) (equal weight 
at ti and t<z), then it is found that no does not contribute at all 
to the determination of fi, and the standard error is still given 
by Eq. (^|). If, on the other hand, the two observation epochs 
are not equivalent, then some improvement can be expected by 
introducing the position measurements. 



An important special case is when there is just a posi- 
tion (no proper motion) determined at one of the epochs, say 
t\. This is however equivalent to the two independent proper- 
motion determinations no at to = (ti + t 2 )/2, and ni at t 2 , 
separated by t 2 — to = T/2. Applying Eq. (Q) on this case 
yields 



e(v r 



2A_ 

TlTfl 



rp2 



1/2 



(7) 



This formula is applicable on the combination of a recent posi- 
tion and proper-motion measurement (e.g. by Hipparcos) with 
a position derived from old photographic plates (e.g. the As- 
trographic Catalogue). Taking t\ ~ 1907, e(</>i) — 200 mas 
as representative for the Astrographic Catalogue, and t 2 = 
1991.25, e{4>2) — 1 mas, e(/i 2 ) ~ 1 mas yr _1 for Hippar- 
cos, we find e(v r ) ~ (60 kms -1 arcsec 2 yr _1 )/(7r/i). With 
such data, moderate accuracies of a few tens of km s _1 can be 



reached for several stars (Sect. 5.3). 



5.2. Effects of gravitational perturbations 

The perspective-acceleration method depends critically on the 
assumption that the star moves with uniform space motion rel- 
ative the observer. The presence of a real acceleration of their 
relative motions, caused by gravitational action of other bod- 
ies, would bias the calculated astrometric radial velocity by 
—gt/2fi, where g t is the tangential component of the rela- 
tive acceleration. The acceleration towards the Galactic centre 
caused by the smoothed Galactic potential in the vicinity of the 
Sun is g ~ 2 x 10 -13 km s~ 2 . For a hypothetical observer 
near the Sun but unaffected by this acceleration, the maximum 
bias would be 0.06 km s _1 for Barnard's star, and 0.17 km s~ 2 
for Proxima. However, since real observations are made rela- 
tive the solar-system barycentre, which itself is accelerated in 
the Galactic gravitational field, the observed (differential) ef- 
fect will be very much smaller. 

In the case of Proxima the acceleration towards a Cen AB 
is of a similar magnitude as the Galactic acceleration. For 
the several orbital binaries in Table ^| the curvature of the or- 
bit is much greater than the perspective acceleration. Applica- 
tion of this method will therefore require careful correction for 
all known perturbations: the possible presence of long-period 
companions may introduce a considerable uncertainty. 

Among other effects which may have to be considered are 
light- time effects which, to first order in c _1 , may require a 
correction of — /2c on the right-hand side of Eq. (Q), where 
Vt is the tangential velocity. For typical high-velocity (Popula- 
tion II) stars the correction is 0.1-0.2 km s _1 . At this accuracy 
level, the precise definition of the radial-velocit y con cept itself 
requires careful consideration (Lindegren et al. 1999). 



5.3. Results from observed proper-motion changes 

Past deter minati ons of the perspective accel eratio n, e.g. by van 
de Kamp ( 1981 ) and Gatewood & Russell ( 1974 ), were based 
on photographic observations collected over several decades, 
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Table 3. Astrometric radial velocities, obtained by combining positions and proper motions from Hipparcos (epoch 1991.25) 
with old position measurements from the Astrographic Catalogue. Spectroscopic radial velocities are also given (from Turon et 
al. |1998| ). Data for the binary 61 Cyg = HIP 104214+104217 = AC 1382645+1382649 refer to its mass centre, assuming a mass 
ratio of 0.90. Two solutions are given: the first based on Hipparcos and AC data alone; the second (marked *) includes Bessel's 
visual measurements from 1838. 



HIP 



AC 2000 



Radial velocity, v r [km s 



Remark 





No. 


Epoch 


Astrom. 


Spectr. 




439 


3152964 


1912.956 


+ * .u 


zt 


9Q 7 


+22.9 




1475 


1406215 


1898.435 


— 40.8 


± 


36.1 


+ 12.0 


GX And 


5336 


1721511 


1913.868 


— 106.4 


± 


82.7 


-98.0 


fi Cas 


15510 


3488626 


1901.018 


+89.6 


± 


59.1 


+86.8 


82Eri 


19849 


2125614 


1892.970 


-73.2 


± 


32.8 


-42.6 


40Eri 


24186 


3505363 


1899.058 


+249.4 


± 


13.7 


+245.5 


Kapteyn's star 


36208 


282902 


1908.859 


-33.7 


± 


38.4 


+18.7 


Luyten's star 


54035 


1340883 


1930.895 


+22.0 


± 


36.1 


-85.6 




54211 


1463341 


1895.620 


+58.6 


± 


30.5 


+67.5 




57367 


4195112 


1924.492 


+43.1 


± 


106.4 






57939 


1342199 


1930.260 


-139.5 


± 


86.1 


-98.0 


Groombridge 1830 


87937 


146626 


1905.979 


-101.9 


± 


6.5 


-109.7 


Barnard's star 


104214/217 


1382645/649 


1921.699 


-12.2 


± 


34.5 


-64.5 


61 Cyg 


104214/217 


1382645/649 


1921.699 


-68.0 


± 


11.1 * 






105090 


3462277 


1905.316 


-25.0 


± 


40.9 


+23.6 




108870 


4384302 


1901.189 


-64.1 


± 


23.6 


-40.0 


elnd 


114046 


3355101 


1913.368 


+54.0 


± 


19.7 


+9.5 





in which the motion of the target star was measured relative 
to several background (reference) stars. One difficulty with the 
method has been that the positions and motions of the reference 
stars are themselves not accurately known, and that small errors 
in the reference data could cause a spurious acceleration of the 
target star (van de Kamp 1935a). 



The Hipparcos Catalogue (ESA 1997 ) established a very 
accurate and homogeneous positional reference frame over the 
whole sky. Using the proper motions, this reference frame can 
be extrapolated backwards in time. It is then possible to re- 
reduce measurements of old photographic plates, and express 
even century-old stellar positions in the same reference frame 
as modern observations. This should greatly facilitate the deter- 
mination of effects such as the perspective acceleration, which 
are sensitive to systematic errors in the reference frame. 

As part of the Carte du Ciel project begun more than a 
century ago, an astrographic programme to measure the posi- 
tions of all stars down to the 11th magnitude was carried out 
and published as the Astrographic Catalogue, AC (see Eich- 



horn 1974 for a description). After transfer to electronic media, 



the position measurements have been reduced to the Hippar- 



cos reference frame (Nesterov et al. 1998, Urban et al. 1998) 



The result is a positional catalogue of more than 4 million stars 
with a mean epoch around 1907 and a typical accuracy of about 
200 mas. We have used the version known as AC 2000 (Urban 



et al. 1998), available on CD-ROM from the US Naval Obser- 



vatory, to examine the old positions of all the stars with HIP 
identifiers in Table ||. 



For the stars in Table || we successfully matched the AC 
positions with the positions extrapolated backwards from the 
Hipparcos Catalogue and hence could calculate the astrometric 
radial velocities. Other potential targets in Table || were either 
outside the magnitude range of AC 2000 (e.g. a Cen and Prox- 
ima) or lacked an accurate proper motion from Hipparcos (e.g. 
van Maanen 2 and HIP 91768+91772). 

The basic procedure was as follows. The rigorous epoch 
transformation algorithm described in Sect. 1 .5.5 of the Hippar- 
cos Catalogue, Vol. 1, was used to propagate the Hipparcos po- 
sition and its covariance matrix to the AC 2000 epoch relevant 
for each star. This extrapolated position was compared with the 
actual measured position in AC 2000, assuming a standard er- 
ror of 200 mas in each coordinate for the latter. A \ 2 goodness- 
of-fit was then calculated from the position difference and the 
combined covariance of the extrapolated and measured posi- 
tions. The epoch transformation algorithm requires that the ra- 
dial velocity is known. The radial velocity was therefore varied 
until the x 2 attained its minimum value. The ±ler confidence 
interval given in the table was obtained by modifying the ra- 
dial velocity until the \ 2 had increased by one unit above the 
minimum. 

For some of the stars, data had to be corrected for duplicity 
or known orbital motion. The solutions for the resolved binary 
61 Cyg (HIP 104214+104217) refer to the mass centre, assum- 
ing a mass ratio of M-q/M^ — 0.90, as estimated by means of 
standard isochrones from the absolute magnitudes and colour 
indices of the components (Soderhjelm, private communica- 
tion). For the astrometric binary ji Cas (HIP 5336) the Hippar- 
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Fig. 2. Comparison of the astrometric radial velocities in Ta- 
ble U with their spectroscopic counterparts. The straight line is 
the expected relationship w r (astrom.) ~ w r (spectr.). The three 
most accurate determinations are indicated. 



cos data explicitly refer to the mass centre using the orbit by 



Heintz & Cantor (1994); the same orbit was used to correct the 



AC position of the primary to the mass centre. No correction 
for orbital motion was used for GX And and 40 Eri. 

Table || gives two solutions for 61 Cyg. The first solution 
was obtained as described above, using only the Hipparcos data 
plus the AC positions for the two components. The second so- 
lution, marked with an asterisk in the tab le, wa s derived by 
including also the observations by Bessel ( 1839 ) from his pi- 
oneering determination of the star's parallax. Bessel measured 
the angular distances from the geometrical centre (half-way be- 
tween the components) of 61 Cyg to two reference stars, called 
a and b in his paper. After elimination of aberration, proper mo- 
tion and parallax, he found the distances 461.6171 ± 0.015 arc- 
sec and 706.2791 ± 0.017 arcsec for the beginning of year 
1838 (B1838.0 = J1838.0022). The uncertainties are our es- 
timates (standard errors) based on the scatter of the residu- 
als in Bessel's solution 'IT. We identified the reference stars 



in AC 2000 and in the Tycho Catalogue (ESA |1997[ ) as a 
= AC 1382543 = TYC 3168 708 1 and b = AC 1382712 = 
TYC 3168 1106 1. Extrapolating the positions from these cat- 
alogues back to B1838 allowed us to compute the position of 
the geometrical centre of 61 Cyg in the Hipparcos/Tycho refer- 
ence frame. This could then be transformed to the position of 
the mass centre, using Bessel's own measurement of the sepa- 
ration and position angle in 61 Cyg and the previously assumed 
mass ratio. Actually, all the available data were combined into 
a x 2 goodness-of-fit measure and the radial velocity was varied 
in order to find the minimum and the ± 1 a confidence interval. 
This gave v r = —68.0 ± 11.1 km s _1 . 



Table ||] also gives the spectroscopic radial velocities when 
available in the literature. A comparison between the astromet- 
ric and spectroscopic radial velocities is made in Fig. |2|. Given 
the stated confidence intervals, the agreement is in all cases 
rather satisfactory. The exercise demonstrates the basic feasi- 
bility of this method, but also hints at some of the difficulties 
in applying it to non-single stars. 

6. Radial velocity from changing angular extent 
(moving-cluster method) 

The moving-cluster method is based on the assumption that the 
stars in a cluster move through space essentially with a com- 
mon velocity vector. The radial-velocity component makes the 
cluster appear to contract or expand due to its changing dis- 
tance (Fig. [[]:). The relative rate of apparent contraction equals 
the relative rate of change in distance to the cluster. This can 
be converted to a linear velocity (in km s _1 ) if the distance 
to the cluster is known, e.g. from trigonometric parallaxes. In 
practice, the method amounts to determining the space veloc- 
ity of the cluster, i.e. the convergent point and the speed of 
motion, through a combination of proper motion and parallax 
data. Once the space velocity is known, the radial velocity for 
any member star may be calculated by projecting the velocity 
vector onto the line of sight. 

The method can be regarded as an in versio n of the classi- 
cal procedure (e.g. Binney & Merrifield 1998 1 by which the 
distances to the stars in a moving cluster are derived from the 
proper motions and (spectroscopic) radial velocities: if instead 
the distances are known, the radial velocities follow. The first 
application of the classical moving -clust er method for distance 
determination was by Klinkerfues ( 1873 ), in a study of the Ursa 
Major system. The possibility to check spectroscopic radial 
velocities against astrometric data was recognised by Klinker- 
fues, but could not then be applied to the Ursa Major cluster due 
to the lack of reliable trigonometric parallaxes. This changed 
with Hertzsprung's (1909) discovery that Sirius probably be- 
longs to the Ursa Major moving group. The relatively large and 
well-determined parallax of Sirius, combined with its consid- 
erable angular distance from the cluster apex, could lead to a 
meaningful estimate for the cluster velocity and hence for the 



radial velocities. Rasmuson (1921) and Smart (1939) appear to 



have been among the first who actually made this computation, 
although mainly as a means of verifying the cluster method 



for distance determination. Later studies by Petrie (1949) and 



Petrie & Moyls (1953) reached formal errors in the astromet- 
ric radial velocities below 1 km s _1 . The last paper concluded 
"There does not appear to be much likelihood of improving the 
present results until a substantial improvement in the accuracy 
of the trigonometric parallaxes becomes possible." 

One of the purposes of the Petrie & Moyls study was to de- 
rive the astrometric radial velocities of spectral type A in order 
to check the Victoria system of spectroscopic velocities. The 



method was also applied to the Hyades (Petrie 1963) but only 
with an uncertainty of a few km s — 1 , Given the expected future 
availability of more accurate proper motions and trigonomet- 
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ric parallaxes, Petrie ( 1962| ) envisaged that one or two mov- 
ing clusters could eventually be used as primary radial-velocity 
standards for early-type spectra. 

Such astrometric data are now in fact available. In Sect. |6.l| 
we derive a rough estimate of the accuracy of the method and 
survey nearby clusters and associations in order to find promis- 
ing targets for its application. An important consideration is to 
what extent systematic velocity patterns in the cluster, in partic- 
ular cluster expansi on, w ill limit the achievable accuracy. This 
is discussed in Sect. 6.2 and Appendix A. In Sect. 5.3 we briefly 
consider the improvement in the distance estimates for individ- 
ual stars resulting from the moving-cluster method. 

The present discussion of the moving-cluster method is 
only intended to highlight its theoretical potential and limita- 
tions. Its actual application requires a more rigorous formula- 
tion, which is developed in a second paper. 

6.1. Potential accuracy 

The accuracy of the astrometric radial velocity potentially 
achievable by the moving-cluster method can be estimated as 
follows. Let b be the (mean) distance to the cluster and con- 
sider a star at angular distance p from the centre of the cluster, 
as seen from the Sun. The projected linear distance of the star 
from the centre of the cluster is b sin p ~ bp, provided the angu- 
lar extent of the cluster is not very large. As the cluster moves 
through space, its linear dimensions remain constant, so that 
pb + pb = 0. Putting p = p (the proper motion relative to the 
cluster centre), b — v r , and b = A /it, gives v r = —Ap/ (pir). 
Now suppose that the parallaxes and proper motions of n clus- 
ter stars are measured, each with uncertainties of e(7r) and e(p). 
Standard error propagation formulae give the expected accu- 
racy in v r as 



e(v r ) ~ A ■ 



<p) 



VrPm 



Ae(p) 



1/2 



(8) 



where p rms is in radians; A is the astronomical unit (Sect. g). 
The expression within the square brackets derives from the un- 
certainty in the mean cluster distance, by which the derived 
radial velocity scales. For the type of (space) astrometry data 
considered here (Case A and B), e{ir) / e(p) is on the order of a 
few years (for Hipparcos the mean ratio is ~ 1.2 yr). The factor 
in brackets can then be neglected except for the most extended 
(and nearby) clusters. 

Under certain circumstances it is not the accuracy of 
proper-motion measurements that defines the ultimate limit on 
e(v r ), but rather internal velocity dispersion among the cluster 
stars. Assuming isotropic dispersion with standard deviation a v 
in each coordinate, one must add <j v n/A quadratically to the 
measurement error e(p) in Eq. (|J). Thus 



e(w r) 



1 



1/2 



Ae(p) 



1/2 



(9) 



is the accuracy achievable for the radial velocity of the clus- 
ter centroid. For the radial velocity of an individual star this 
uncertainty must be increased by the internal dispersion. 

The internal velocity dispersion will dominate the error 
budget for nearby clusters, viz. if tt > Ae(p)/a v . Assuming 
a velocity dispersion of 0.25 km s -1 and a proper-motion ac- 
curacy of 1 mas yr -1 (as for Hipparcos), this will be the case 
for clusters within 50 pc of the Sun. For an observational ac- 
curacy in the 1-10 pas yr -1 range the internal dispersion will 
dominate in practically all Galactic clusters and Eq. (||) can be 
simplified to e(vo r ) ~ cr v /(p lms y/n). In this case the achiev- 
able accuracy becomes independent of the astrometric one. 

Table ^| lists some nearby clusters and associations, with es- 
timates of the achievable accuracy in the radial velocity of the 
cluster centroid, assuming current (Hipparcos-type) astromet- 
ric performance (Case A in Sect. ^ as well as future (microarc- 
sec) expectations (Case B). As explained above, increasing the 
astrometric accuracy still further gives practically no improve- 
ment; this is why Case C is not considered in the table. 

The entry 'HIP 98321' refers to the possible association 
identified by Platais et al. (1998) and named after one of its 
members. Of dubious status, it was included as an example of 
the extended, low-density groups that may exist in the general 
stellar field, but are difficult to identify with existing data. 

6.2. Internal velocity fields, including cluster expansion 



Blaauw ( 1964 ) showed that the proper motion pattern for a 
linearly expanding cluster is identical to the apparent conver- 
gence produced by parallel space motions. Astrometric data 
alone therefore cannot distinguish such expansion from a ra- 
dial motion. If such an expansion exists, and is not taken into 
account in estimating the astrometric radial velocity, a bias will 
result, as examined in Appendix A (Eq. A. 7). 

The gravitationally unbound associations are known to 
expand on timesca les c omparable with their nuclear ages 
(de Zeeuw & Brand 1985 ). But also for a gravitationally bound 
open cluster some expansion can be expected a s a res ult of the 
dyn amica l evolution of the cluster (see Mathieu 1985 and Wie- 
len 1988 for an introduction to this complex issue). In either 
case the inverse age of the cluster or association may be taken 
as a rough u pper limit on the cluster's relative expansion rate k 
[yr -1 ]. Eq. (A. 7) then gives 



I W«br)| £ 0.9543 x 



fro [pc] 
age [Myr] 



km s 



(10) 



for the bias of a star near the centre of the cluster. (For an 
expanding cluster, (5 cxp is always negative.) Resulting values, 
in the last column of Table [|, are adequately small for a few 
nearby, relatively old clusters. In other cases the potential bias 
is very large and will certainly limit the applicability of the 
method. The OB associations are particularly troublesome, 
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Table 4. Potential accuracy in astrometric radial velocities of nearby star clusters and associations, using the moving-cluster 
method. The first five data columns contain general cluster data taken from the literature (see below): n = number of member 
stars; (photometric) ages; p lms = rms angular radius; 60 = mean distance from the Sun; v r , Vt = approximate radial and 
tangential velocities. The columns headed e(uo r ) give the potential accuracies according to Eq. assuming e(7r) = 1 mas 
and e(/i) = 1 mas yr _1 (Case A), or e(iv) — 1 /ias and e(fi) = 1 /ias yr _1 (Case B). The final column gives the bias in the 
astrometric radial velocity estimate (Eq. [!(]) that would result from neglecting an isotropic expansion of the cluster, assuming that 
the relative expansion rate equals the inverse age of the cluster. A dash means that no data were available. The cluster data are 
mostly from Lynga (1987a 1987b) [assuming p lms to equal the cluster core diameter D in that catalogue], with supplementary 



data fr om de Zeeuw et al. ( 1999|), de Zeeuw & Brand ( |1985[ ), Mermilliod et al. ( |1996| ), Perryman et al. ( |1998[ ), Smith & Messer 



([1983]), Soderblom & Mayor (|1993|), and Stryker & Hrivnak (|1984j). Tangential velocities and the angular radii of the large r 
clusters and of all the associations were calculated from data for identifie d mem bers in the Hippar cos C atalogue (ESA 1997 ). 
The memberships for the associations were adopted from de Zeeuw et al. ( 1999 1 and Platais et al. ( 1998 ). Ursa Major refers to 
the rather exten ded gr oup (UMaG) identified by Soderblom & Mayor ( 1993 ); 'HIP 98321' to the possible association identified 
by Platais et al. ( 1998 ). An internal velocity dispersion a v = 0.25 km s _1 was assumed for all clusters and associations, although 
a higher value is more likely for associations and the halos of some clusters. 



Name 



IAU 
designation 



Age 
[Myr] 



Prms 

[arcmin] 



60 

[pc] 



V r Vt 

[kms -1 ] 



e(v 0r ) [km s 1 ] 
(A) (B) 



5exp(V0r) 

[kms -1 ] 



Cassiopeia-Taurus 




83 


25 


1800 


190 


+6 


21 


0.24 


0.06 


-7.3 


Upper Centaurus Lupus 




221 


13 


670 


140 


+5 


21 


0.25 


0.09 


-10 


Ursa Major 




40 


300 


4300 


25 


-11 


5 


0.11 


0.10 


-0.08 


Lower Centaurus Crux 




180 


10 


560 


118 


+ 12 


19 


0.30 


0.12 


-11 


Hyades 


t, 1)424+ 15/ 


1 on 
380 


0/5 


301) 


40 


1 A 1 

+43 


Id 


n 1 n 


A 1 A 

V. 14 


n (yi 
— U.U / 


Perseus OB 3 (a Per) 


C 0318+484 


186 


50 


350 


180 


-1 


29 


0.64 


0.18 


-3.4 


'HIP 98321' 




59 


60 


740 


300 


-17 


4 


1.1 


0.19 


-4.8 


Upper Scorpius 




120 


5 


325 


145 


-5 


18 


0.71 


0.24 


-28 


Lacerta OB 1 




96 


16 


350 


370 


-13 


8 


1.8 


0.26 


-22 


Collinder 121 




103 


5 


290 


540 


+26 


15 


3.3 


0.32 


-103 


Collinder 70 


C 0533-011 


345 




140 


430 






2.7 


0.33 




Cepheus OB2 




71 


5 


320 


615 


-21 


12 


4.1 


0.34 


-120 


Vela OB2 




93 


20 


260 


415 


+ 18 


20 


2.8 


0.36 


-20 


Perseus OB2 




41 


7 


340 


300 


+20 


14 


2.5 


0.43 


-46 


Pleiades 


C 0344+239 


277 


130 


120 


125 


+7 


29 


1.1 


0.43 


-0.92 


Coma Berenices 


C 1222+263 


273 


460 


120 


87 





9 


0.84 


0.43 


-0.18 


NGC 3532 


C 1104-584 


677 


290 


50 


480 


+7 


27 


6.0 


0.66 


-1.6 


Praesepe 


C 0837+201 


161 


830 


70 


160 


+33 


26 


3.1 


0.98 


-0.18 


NGC 2477 


C 0750-384 


1911 


1260 


20 


1150 


+7 


17 


21 


0.98 


-0.87 


IC 4756 


C 1836+054 


466 


830 


39 


390 


-18 


5 


7.6 


1.0 


-0.45 


IC 4725 


C 1828-192 


601 


41 


29 


710 


+3 


20 


16 


1.2 


-17 


Trumpler 10 




23 


15 


150 


370 


+21 


26 


8.6 


1.2 


-24 


Cepheus OB 6 




20 


50 


150 


270 


-20 


22 


6.8 


1.3 


-5.2 


NGC 752 


C 0154+374 


77 


3300 


75 


360 


-3 


20 


9.0 


1.3 


-0.10 


NGC 6618 


C 1817-162 


660 




25 


1500 




26 


38 


1.3 




NGC 2451 


C 0743-378 


153 


41 


50 


315 


+27 


27 


8.4 


1.4 


-7.3 


NGC 7789 


C 2354+564 


583 


1780 


25 


1800 


-54 


27 


49 


1.4 


-1.0 


NGC 2099 


C 0549+325 


1842 


200 


14 


1300 


+8 


45 


35 


1.4 


-6.4 


NGC 6475 


C 1750-348 


54 


130 


80 


240 


-12 


7 


6.8 


1.5 


-1.8 


NGC 2264 


C 0638+099 


222 


10 


39 


800 


+22 


17 


23 


1.5 


-76 


Stock 2 


C 0211+590 


166 


170 


45 


300 


+2 


30 


8.6 


1.5 


-1.7 


IC 2602 


C 1041-641 


33 


29 


100 


150 


+22 


14 


4.6 


1.5 


-5.0 



not only because they are young objects (implying large val- 
ues of k), but also because they sometimes appear to expand 
significantly faster than their photometric ages would suggest 
(de Zeeuw & Brand |1985|). 



However, it should be remembered that the ultimate limi- 
tation set by cluster expansion depends on how accurately the 
expansion rate k can be estimated by some independent means. 
For instance, if k can somehow be estimated to within 10 per 
cent of its value, then the residual biases would still be on the 
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sub-km s" 1 level for most of the objects in Table [I| Numerical 
simulation of the dynamical evolution of clusters might in prin- 
ciple provide such estimates of k, as could the spectroscopic ra- 
dial velocities as function of distance. The use of spectroscopic 
data would not necessarily defeat the purpose of the method, 
i.e. to determine absolute radial velocities, since the expansion 
is revealed already by relative measurements. 

6.3. Distances to the individual stars 

In a rigorous estimation of the space motion of a moving clus- 
ter, such as will be presented in a second paper, the distances to 
the individual member stars of the cluster appear as parameters 
to be estimated. A by-product of the method is therefore that 
the individual distances are improved, sometimes considerably, 
comp ared w ith the o rigina l trigonometric distances (Dravins 
et al. 1997; Madsen 1999). The improvement results from a 
combination of the trigonometric parallax 7r tr i g with the kine- 
matic (secular) parallax 7Tki n = A[i/v t derived from the star's 
proper motion /i and tangential velocity v t , the latter obtained 
from the estimated space velocity vector of the cluster. The ac- 
curacy of the combined parallax estimate tt can be estimated 
from e(7r)~ 2 = e(7r tr ig) _2 + e(7r kin )~ 2 . In calculating e(irun) 
we need to take into account the observational uncertainty in fj, 
and the uncertainty in v t from the internal velocity dispersion. 
The result is 



6(7?) 



eKrig) 2 + 



1-1/2 



e( M ) 2 + (na v /Ay 



(ID 



From the data in Table |4| we find that, in Case A, the moving- 
cluster method will be useful to resolve the depth structures of 
the Hyades cluster and of the associations Cassiopeia-Taurus, 
Upper Centaurus Lupus, Lower Centaurus Crux, Perseus OB3 
and Upper Scorpius. In Case B, all the clusters and associations 
are resolved by the trigonometric parallaxes, so the kinematic 
parallaxes will bring virtually no improvement. 

Calculation of kinematic distances to stars in moving clus- 
ters is of course a classical procedure (e.g., Klinkerfues 1873 



and van Bueren 1952); what is new in our treatment is that such 



distances are derived without recourse to spectroscopic data. 



7. Further astrometric methods 

With improved astrometric data, further methods for radial- 
velocity determinations may become feasible. The moving- 
cluster method could in principle be applied to any geomet- 
rical configuration of a fixed linear size. To reach an accuracy 
of 1 km s _1 in the astrometric radial velocity of an object at 
10 pc distance requires a dimensional 'stability' of the order of 
10~ 7 yr _1 ; at a distance of 1 kpc the requirement is 1CP 9 yr _1 . 
Since these numbers are greater than or comparable with the in- 
verse dynamical timescales of many types of galactic objects, 
there is at least a theoretical chance that the method could work, 
given a sufficient measurement accuracy. We consider briefly 
two such possibilities. 



7.1. Binary stars 

According to the previous argument it would be possible to ig- 
nore the relative motion of the components in a binary if the 
period is longer than some 10 7 yr. This implies an linear sepa- 
ration of at least some 50 000 astronomical units, or a few de- 
grees on the sky at 10 pc distance. In principle, then, this case 
is equivalent to a moving cluster with n — 2 stars. 

In the opposite case of a (relatively) short-period binary, the 
radial velocity might be obtained from apparent changes of the 
orbit. The projected orbit will not be closed, but form a spiral on 
the sky: slightly diverging if the stars are approaching, slightly 
converging if they recede. For a system at a distance of 10 pc, 
say, with a component separation of 10 astronomical units, a 
radial velocity of 100 km s _1 will change the apparent orbital 
radius of 1 arcsec by 10 /xas per year. The relative positions 
would need to be measured during at least a significant fraction 
of an orbital period, or some 20 years in our example, resulting 
in an accumulated change by about 0.2 mas. 

Since only relative position measurements between the 
same stars are required, the observational challenges are not 
as severe as in some other cases. For a binary with a few arcsec 
separation, the isoplanatic properties of the atmosphere imply 
that the cross-correlation distance between the speckle images 
of both stars should be stable to better than one mas. Averag- 
ing very many exposures should reduce the errors into the /ias 
range, with pra ctical limits possibly set by differential refrac- 
tion (McAlister [l996| ). 



7.2. Globular clusters 

The moving-cluster method (Sect. ^) could in principle be 
applied also to globular clusters. Globular clusters differ 
markedly from open clusters in that (potentially) many more 
stars could be measured, and through a m uch la rger velocity 
dispersion (~ 5 km s^ 1 ; Peterson & King 1975 ). The higher 
number of stars partly compensates the larger dispersion. How- 
ever, all globular clusters are rather di stan t, making their angu- 
lar radii small. As discussed in Sect. 5.1 the approximate for- 
mula e(for) — cr v / (PrmsV^) applies in the case when the in- 
ternal motions are well resolved. Taking p rms ~ 10-20 arcmin 
as representative for the angular radii of comparatively nearby 
globular clusters, we find that averaging over some 3 x 10 4 to 
10 5 member stars is needed to reach a radial-velocity accuracy 
comparable with a v , i.e. several km s _1 . Furthermore, in view 
of the discussion in Sect. 6.2, it is not unlikely that the complex 
kinematic structures of these objects (e.g. Norris et al 



1997) 



would bias the results. Thus, globular clusters remain difficult 
targets for astrometric radial-velocity determination. 



8. Conclusions 

The theoretical possibility to use astrometric data (parallax and 
proper motion) to deduce the radial motions of stars has long 
been recognised. With the highly accurate (sub-milliarcsec) as- 
trometry already available or foreseen in planned space mis- 
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sions, such radial-velocity determinations are now also a prac- 
tical possibility. This will have implications for the future def- 
inition of radial-velocity standards, as the range of geometri- 
cally determined accurate radial velocities, hitherto limited to 
the solar system and to solar-type spectra, is extended to many 
other stellar types represented in the solar neighbourhood. 

We have identified and analysed three main methods to 
determine astrometric radial velocities. The first method, us- 
ing the changing annual parallax, is the intuitively most obvi- 
ous one, but requires data of an accuracy beyond current tech- 
niques. It is nevertheless potentially interesting in view of fu- 
ture space missions or long-term observations from the ground. 

The second method, using the changing proper motion or 
perspective acceleration of stars, has a long history, and was 
previously applied to a few objects, albeit with modest preci- 
sion in the resulting radial velocity. New results for a greater 
number of stars, obtained by combining old and modern data, 
were given in Table |] and Fig. 0, thus proving the concept. 
However, to realise the full potential of the method again re- 
quires the accuracies of future astrometry projects. 

In both these methods the uncertainty in the astrometric 
radial velocity increases, statistically, with distance squared. 
They are therefore in practice limited to relatively few stars 
very close to the Sun and, in the second method, to stars with a 
large tangential velocity. In the general case, the two methods 
could actually be combined to yield a somewhat higher accu- 
racy, but at least for the stars considered in Tables [ljand ^Jthis 
would only bring a marginal improvement. 

The third method, using the changing angular extent of a 
moving cluster or association, is an inversion of the classical 
moving-cluster method, by which the distance to the cluster 
was derived from its radial velocity and convergence point. If 
the distance is known from trigonometric parallaxes, one can 
instead calculate the radial velocities. It appears to be the only 
method by which astrophysically interesting accuracies can be 
obtained with existing astrometric data. In future papers we 
will develop and exploit this possibility in full, using data from 
the Hipparcos mission. A by-product of the method is that the 
distance estimates to individual cluster stars may be signifi- 
cantly improved compared with the parallax measurements. 

One would perhaps expect the moving-cluster method to 
become extremely powerful with the much more accurate data 
expected from future astrometry projects. Unfortunately, this 
is not really the case, as internal velocities (both random and 
systematic) become a limiting factor as soon as they are re- 
solved by the proper-motion measurements. Nevertheless, even 
the limited number of clusters within reach of such determi- 
nations contain a great many stars spanning a wide range in 
spectral type and luminosity. 
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Appendix A: Effects of internal velocities on the 
moving-cluster velocity estimate 

In this Appendix we examine how sensitive the moving-cluster 
method is to systematic velocity patterns in the cluster, and to 
what extent such patterns can be determined independently of 
the astrometric radial velocity. For this purpose we may ignore 
the random motions as well as the observational errors and we 
consider only a linear (first-order) velocity field. 

Let 60 be the position of the cluster centroid relative to the 
Sun and s = b — bo the position of a member star relative to 
the centroid. The space velocity of the star is v = Vq + 77, 
where r\ is the peculiar velocity. The velocity field is described 
by the tensor T such that 77 = Ts. In Cartesian coordinates 
the components of this tensor are simply the partial deriva- 
tives dv a /dbj3 for a,{3 = x, y, z. These nine numbers together 
describe the three components of a rigid-body rotation, three 
components of an anisotropic expansion or contraction, and 
three components of linear shear. 

It is intuitively clear that certain components of the linear 
velocity field, such as rotation about the line of sight, can be 
determined purely from the astrometric data. If the correspond- 
ing components of T are included as parameters in the cluster 
model, they can be estimated and will not produce a systematic 
error (bias) in the astrometric radial velocity derived from the 
model fitting. Such components of T are 'observable' and in 
principle not harmful to the method. Let us now examine more 
generally the extent to which T is observable by astrometry. 

Suppose there exists a non-zero tensor T such that the ve- 
locity fields Vq + Ts and Vq + Av produce identical obser- 
vations for some vector Av. Since the cluster velocity Vo is 
already a parameter of the model, the observational effects of 
the velocity field T could then entirely be absorbed by adjust- 
ing Vq. In this case T would be a non-observable component of 
the general velocity field. Moreover, if there exists such a com- 
ponent in the real velocities, then the estimated cluster velocity 
will have a bias equal to Av. 

We now need to calculate the effect of the arbitrary field T 
on the observables. Since the parallaxes are not affected, only 
the proper motion vector 

r = (I - rr')vn/A (A.l) 

needs to be considered. In this equation r is the unit vector from 
the Sun towards the star, f is the rate of change of that direction, 
and / — rr' is the tensor representing projection perpendicular 
to r [I is the unit tensor; thus (J — rr')x = x — rr'x is the 
tangential component of the general vector x]. With b = \b\ = 
A/trwe can write s = rb— bo. T is non-observable if the space 
velocities Vq + Ts and Vq + Av produce identical tangential 
velocities for every star, i.e. if 

(I - rr')[v + T(rb - bo)] = (I - rr')(«o + Aw) (A.2) 
for all directions r and distances b. This is equivalently written 
(I - rr')Trb -(I- rr')(Tb + Av) = . (A.3) 
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In order that this should be satisfied for all b, it is necessary that 
(I-rr')Tr = and (I-rr')(Tb + Av) = are separately 
satisfied for all unit vectors r. The latter equation implies that 



Av = -Th 



o ■ 



(A.4) 



The former equation can be written Tr = rr'Tr, which shows 
that r is an eigenvector of T (with eigenvalue r'Tr). But the 
only tensor for which every unit vector is an eigenvector is the 
isotropic tensor, T = Ik for the arbitrary scalar k. It follows 
that the only non-observable component of T is of the form 
Ik, parametrised by the single scalar k, and that consequently 
eight linearly independent components of T can, in princi- 
ple, be determined from the astrometric observations. The non- 
observable field T — Ik describes a uniform isotropic expan- 
sion (k > 0) or contraction (k < 0) of the cluster with respect 
to its centroid. These effects are observationally equivalent to 
an approach or recession of the cluster, i.e. to a different value 
of its radial velocity. 

Av is the bias for the centroid velocity. For any given star, 
the bias vector Sv is the difference between the derived (appar- 
ent) space velocity vector Vq + Av and the true vector Vq+Ts. 



Using s = b — bo and Eq. ( A.4 ) we find 



Sv = Tb . 



(A.5) 



The resulting bias in the astrometric radial velocity is 

5{v r ) = -r'Tb. (A.6) 

Isotropic expansion (T = Ik), in particular, gives the bias 

ScxpM = -bit. (A.7) 

For a uniformly expanding cluster k^ 1 equals the expansion 
age, i.e. the time elapsed since all the stars were confined to a 
very small region of space. 
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